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The main goal of this work is to determine which entire functions preserve nonnegativity of 
matrices of a fixed order n — i.e., to characterize entire functions / with the property that f{A) is 
entrywise nonnegative for every entrywise nonnegative matrix A of size nxn. Towards this goal, 
we present a complete characterization of functions preserving nonnegativity of (block) upper- 
triangular matrices and those preserving nonnegativity of circulant matrices. We also derive 
necessary conditions and sufficient conditions for entire functions that preserve nonnegativity 
of symmetric matrices. We also show that some of these latter conditions characterize the even 
or odd functions that preserve nonnegativity of symmetric matrices. 

1 Motivation 

The purpose of this paper is to investigate which entire functions preserve nonnegativity of matrices 
of a fixed order. More specifically, we consider several classes of structured matrices whose struc- 
ture is preserved by entire functions and characterize those entire functions / with the property 
that f{A) is entrywise nonnegative for each entrywise nonnegative matrix A of size nxn. The 
characterizations that we obtain might be of independent interest in matrix theory and other areas 
of mathematics. One of our own motivations behind our investigation is its relevance to the inverse 
eigenvalue problem for nonnegative matrices. 

The long-standing inverse eigenvalue problem for nonnegative matrices is the problem of de- 
termining, given an n-tuple (multiset) A of complex numbers, whether there exists an entrywise 
nonnegative matrix A whose spectrum cr{A) is A. The literature on the subject is vast and we 
make no attempt to review it. The interested reader is referred to books [25] and [1], expository 
papers [3], [9], [18], [19] and references therein, as well as to some recent papers [23], [1], [21], [30] . 



The necessary conditions for a given n-tuple to be realizable as the spectrum of a nonnegative 
matrix known so far for arbitrary values of n can be divided into three groups: conditions for 
nonnegativity of moments, Johnson-Loewy-London inequalities, and Newton's inequalities. 
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Given an n-tuple A, its moments are defined as follows: 

s^(A): = ^A"^, me IN. 

AeA 

If A = cr{A) for some nonnegative matrix A, then Sm(A) is nothing but the trace tr{A"^), and 
therefore must be nonnegative. Another basic condition follows from the Perron-Frobenius the- 
ory [27], [H]: the largest absolute value max;veA 1^1 must be the Perron eigenvalue of a realizing 
matrix A and therefore must itself be in A. Finally, the multiset A must be closed under complex 
conjugation, being the spectrum of a real matrix A. Interestingly, the last two conditions are in 
fact not independent conditions, but follow from the nonnegativity of moments, as was shown by 
Friedland in [10]. Thus, there turns out to be just one set of basic conditions 

s-m(A) > for m G IN. 

The next set of necessary conditions was discovered independently by Loewy and London in [22j 
and by Johnson in [17] . These conditions relate moments among themselves as follows: 

s^(A) < n™-isfc„(A), A;,mG IN. 

Newton's inequalities were conjectured in [14] and proved for M-matrices in |13j . An M-matrix 
is a matrix of the form rl — A, where ^ is a nonnegative matrix, r > q{A), and where g{A) is the 
spectral radius of A: 

q{A) : = max |A|. 

AGcr(A) 

If M is an M-matrix of order n, then the normalized coefficients Cj{M) of its characteristic poly- 
nomial defined by 

det(A/ - M) =: ^(-l)-'' Cj(M)A"-^' 
i=o 

must satisfy Newton's inequalities 

c]{M) > c,„i(M)c,+i(M), j = 1, . . . ,n - 1. 

Since the coefficients Cj{M) are determined entirely by the spectrum of M, and the latter is obtained 
from the spectrum of a nonnegative matrix A by an appropriate shift, Newton's inequalities form yet 
another set of conditions necessary for an n-tuple to be realizable as the spectrum of a nonnegative 
matrix. The above three sets of conditions — i.e., nonnegativity of moments, Johnson-Loewy- 
London inequalities and Newton's inequalities are all independent of each other but are not sufficient 
for realizability of a given n-tuple (see [13]). 

Quite a few sufficient conditions are also known (see, e.g., [37], [19], [lOj, ^) as well as certain 
techniques for perturbing or combining realizable n-tuples into new realizable n- or m-tuples (where 
m > n) (see, e.g., [M], [33], [H]). Also, necessary and sufficient conditions on an n-tuple to serve as 
the nonzero part of the spectrum of some nonnegative matrix are due to Boyle and Handelman [2j . 

Finally, it follows from the Tarski-Seidenberg theorem [38^ i29j that all realizable n-tuples form 
a semialgebraic set (see also [E]), i.e., for any given n, there exist only finitely many polynomial 
inequalities that are necessary and sufficient for an n-tuple A to be realizable as the spectrum of 
some nonnegative matrix A (this observation was communicated to us by S. Friedland): 
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Indeed, each realizable n-tuple A = (Ai, . . . , A^) is characterized by the condition 

n 

3A>Q : det(A/-^) = JJ(A- Aj). 

i=i 

The last condition is equivalent to each elementary symmetric function cTj(A) being equal to the 
jth coefficient of the characteristic polynomial of A multiplied by (—!)■' — i.e., to the sum of all 
principal minors of A of order j, for j = 1, . . . ,n. Since the set of all nonnegative matrices is a 
semialgebraic set in m? entries of the matrix and since each sum of all principal minors of A of 
order j is a polynomial in the entries of A, the lists of coefficients of characteristic polynomials of 
nonnegative matrices form a semialgebraic set, and hence the n-tuples whose elementary symmetric 
functions match one of those lists also form a semialgebraic set by the Tarski-Seidenberg theorem. 

However, despite so many insights into the subject, and despite the results obtained so far, 
the nonnegative inverse eigenvalue problem remains open. In fact, the problem remains open when 
specialized to several important classes of structured matrices — for instance, the class of entrywise 
nonnegative symmetric matrices. 

Note that the three sets of conditions on an n-tuple A that we discussed above — i.e., non- 
negativity of moments, the Johnson-Loewy-London inequalities and the Newton inequalities — 
are necessary conditions for the realizability of A as the spectrum of a symmetric n x n matrix 
with nonnegative entries (provided, of course, that all the entries of A are now real). A signifi- 
cant fraction of this paper will be devoted to an idea that has relevance to the inverse eigenvalue 
problem for nonnegative symmetric matrices. It is an idea that was first expressed by Loewy and 
London in p2]. When adapted to symmetric matrices, it may be stated as follows: Suppose a 
primary matrix function / is known to map nonnegative symmetric matrices of some fixed order 
n into themselves. Thus f{A) is nonnegative whenever A is. Since f{a{A)) = a{f{A)), both the 
spectrum cr{A) and its image under the map / must then satisfy the aforementioned conditions for 
realizability. This enlarges the class of necessary conditions for the symmetric nonnegative inverse 
eigenvalue problem. Describing this larger class would require knowing exactly what functions / 
preserve nonnegativity of such matrices matrices (of a fixed order). Towards this end, we provide 
a characterization of all the even and odd entire functions that preserve entrywise nonnegativity of 
nonnegative symmetric matrices. 

Along the way, we also obtain complete characterizations of all entire functions that preserve 
nonnegativity of the following classes of structured matrices: 

• Triangular and block-triangular matrices 

• Circulant matrices 

We ought to add here that, for the above classes of structured matrices, our results do not have a 
bearing on the nonnegative inverse eigenvalue problems associated to them. In fact, the solutions 
of the latter problems are quite straightforward. To be precise: an n-tuple A is the spectrum of an 
n X n triangular matrix if and only if all the entries of A are non-negative. As for circulants: the 
eigenvalues of a circulant matrix A are determined by its first row a := [oq ai . . . a„_i] (see [7J), and 
in fact, there is a constant matrix W (i.e., independent of a and A) such that cr{A) = aW. Thus 
the realizable n-tuples in this case are of the form aW, a € IR". Nevertheless, we feel that the 
problem of characterizing the functions that preserve nonnegativity of the above classes of matrices 
can be of interest, independent of the nonnegative inverse eigenvalue problem. 
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2 Outline 



This paper is organized as follows. We make several preliminary observations in Section HI Be- 
fore focusing attention on aspects of the symmetric nonnegative inverse eigenvalue problem, we 
study the structured matrices just discussed. In section [H we characterize the class of functions 
preserving nonnegativity of triangular and block triangular matrices. It turns out that these are 
characterized by nonnegativity conditions on their divided differences over the nonnegative reals. 
Next, in Section [6l we obtain a characterization of functions preserving nonnegativity of circulant 
matrices. This characterization is quite different from that in Section [5] — it involves linear combi- 
nations of function values taken at certain non-real points of C In Section [TJ we obtain a complete 
characterization of the class Tn for small values of n. 

The remainder of the paper is essentially devoted to functions that preserve nonnegativity of 
symmetric matrices. In Section 18. H we review existing results in that direction. In particular, 
we discuss the restriction of [24j, Corollary 3.1] to entire functions, which claims to provide a 
characterization of entire functions that preserve entrywise nonnegativity of symmetric matrices of 
a fixed order. We point out that, while this result is true when restricted to nonnegative definite 
nonnegative symmetric matrices, the condition occurring in that result is not sufficient for an 
entire function to preserve nonnegativity of all symmetric matrices. The techniques leading to [24^ 
Corollary 3.1], however, turn out to be very useful. We use these techniques, along with some new 
ideas, to obtain necessary conditions and sufficient conditions, and characterizations of the even 
and odd entire functions that preserve nonnegativity of symmetric matrices of a fixed order. This 
is the content of Sections 18.21 and 18.31 Because of a gap between the necessary and the sufficient 
conditions, which we also point out in Section 18.21 the results of that section do not provide a 
characterization of all functions preserving nonnegativity of symmetric matrices. We end the paper 
with a list of several open problems in Section [U and suggest various approaches to their solution 
that we have not explored in this paper. 

3 Notation 

We use standard notation IR'"^" for real matrices of size mxn, 1R-|- for nonnegative reals and 
for nonnegative integers, A > > 0) to denote that a matrix A is entrywise nonnegative 
(positive), and cr{A) to denote the spectrum of A. For x G IR, we use [x\ to denote the greatest 
integer that is less than or equal to x. 

4 Preliminaries 

The main goal of the paper is to characterize functions / such that the matrix f{A) is (entrywise) 
nonnegative for any nonnegative matrix A of order n. Since the primary matrix function f{A) is de- 
fined in accordance with values of / and its derivatives on the spectrum of A (see, e.g., |15t Sections 
6.1, 6.2]), we want to avoid functions that are not differentiable at some points in C Therefore we 
restrict ourselves to functions that are analytic everywhere in C, i.e., to entire functions. Thus we 
consider the class 

J'n-={f entire : A G IR"^", ^ > ^ f{A) > 0}. 
Note right away that the classes J-n are ordered by inclusion: 
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Lemma 1 For any n G IN, ^ ^n+i ■ 

Proof. Let A be a nonnegative matrix of order n and let / € Tn+i- Consider the block- 
diagonal matrix B : = diag(A, 0) obtained by adding an extra zero row and column to A. Since 
f{B) = diag(/(A), 0), the matrix f{A) must be nonnegative. Thus / G J^n- D 

Recall that any entire function can be expanded into its Taylor series around any point in C, 
and that the resulting series converges everywhere (see, e.g., [5]). We will mostly focus on Taylor 
series of functions in !Fn centered at the origin. We start with some simple observations regarding 
a few initial Taylor coefficients of such a function. 



Proposition 2 Let f{z) = X^^^q ajz^ be a function in Tn- Then, aj > for j = 0, . 



,n 



1. 



Proof. For n = 1, the statement follows from evaluating / at 0. If n > 1 and / € J^n, then 
evaluate the function / at the matrix 



A:. 
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the entries oq, . . ., fln-i of f{A) must be nonnegative. This finishes the proof. 

Corollary 3 A function f is in Tn for all n (^IN if and only if it has the form f{z) 
with Oj > for all j E . 



Proof. One direction follows from Proposition O The other direction is trivial: if all terms 
in the Taylor expansion of / around the origin are nonnegative, then f{A) combines powers of a 
nonnegative matrix A using nonnegative coefficients, so the resulting matrix is nonnegative. Here 
we make use of the standard fact [151 Theorem 6.2.8] that the matrix power series Yl'jLo ^i^"' 
converges to f{A). □ 

Remark. It must be noted that Proposition [2] cannot be a necessary condition for an entire 
function to belong to J^n- This is easy to see; fix an n E IN and set 



F{x) 



n-l 
j=0 
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where we choose aj > 0, j = 0, . . . , n — 1. Then, there exists an xq > such that F{x) < Vx G 
(xo, cxd). If we set A = rl, for some r € (xq, oo), then A is entrywise nonnegative while the diagonal 
entries of F{A) are negative. Hence, although aj > for j = 0, . . . ,n — 1, F does not preserve 
nonnegativity. 

To conclude this section, we make two more general observations. 

Lemma 4 An entire function f belongs to Tn if and only if it maps positive matrices of order n 
into nonnegative matrices. 

Proof. This is simply due to the continuity of /, since the set of strictly positive matrices is 
dense in the set of all nonnegative matrices of order n. □ 

Lemma 5 For any primary matrix function f , any permutation matrix P and any diagonal matrix 
D with positive diagonal elements, f{A) is nonnegative if and only if f{PDA{PD)~^) is nonnega- 
tive. 

Proof. Note that {P D) f {A){P D)"^ = f{PDA{PD)-^) and that both matrices PD and {POy^ 
are nonnegative. So, f{A) is nonnegative if and only if the matrix f {P D A{P D)"^) is nonnegative. 
□ 

We now analyze three superclasses of our class .Fn • 

• entire functions preserving nonnegativity of upper-triangular matrices; 

• entire functions preserving nonnegativity of circulant matrices; 

• entire functions preserving nonnegativity of symmetric matrices. 

5 Preserving nonnegativity of (block-)triangular matrices 

We first discuss functions preserving nonnegativity of upper- (or lower-) triangular matrices. The 
characterization that we obtain makes use of the notion of divided differences. The divided difference 
(see, e.g., |8]) of a smooth function / at points xi, . . ., Xfc (which can be thought of as an ordered 
sequence xi < ■ ■ ■ < Xfc) is usually defined via the recurrence relation 



and where /[x] :=/(x). Divided differences play a large part in this paper. We shall, however, 
make no attempt to review the results on divided differences that we shall draw upon, especially 
since they are quite readily accessible. The interested reader is referred to [H]. 

Theorem 6 An entire function f preserves nonnegativity of upper-triangular matrices of order n 
if and only if its divided differences of order up to n are nonnegative over IR+, i.e., 




/[xi, . . . , Xfc] > for Xl, . . . , Xfc > 0, k = l,...,n, 



(1) 



or, equivalently, that all derivatives of f of order up to n—1 are nonnegative on 21+. 
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Proof. Sufficiency: Let A=:{aij) be a nonnegative upper-triangular matrix. Suppose a function 
/ satisfies (JI]). By [28], [35] (see also [36]), the elements of the matrix f{A) can be written explicitly 
as 

fiau) i = h 

i>j. 



The divided differences appearing in the sum on the right-hand side are of order not exceeding n; 
hence all the summands, and therefore the sums, are nonnegative. 

Necessity: We proceed by induction on n. If / preserves nonnegativity of upper-triangular 
matrices of order n, it does so also for matrices of order n — 1. Thus, by our inductive hypothesis, 
([l]) holds up to order n — 1. To see that all divided differences of order n are also nonnegative over 
nonnegative reals, consider the matrix A whose first upper diagonal consists of ones, main diagonal 
of n arbitrary nonnegative numbers xi, . . . and all of whose other entries are zero. Then, ([2]) 
shows that f{A)in = f[xi, . . . , x„] and must be nonnegative. 

Finally, since all divided differences of a fixed order k at points in a domain D are nonnegative if 
and only if f^''~^\x) is nonnegative for every point x E D [8J, we see that condition ([T]) is equivalent 
to all derivatives of / of order up to n—1 being nonnegative on ]R-|-. This finishes the proof. □ 

The proofs of ([2]) in [35] and [28] are based on the following observation. 
Result 7 ([35], [28]) A block-triangular matrix of the form 



M 



A B 
a 



a G C\a{A), 



is mapped to the matrix 



fiM) 



f{A) {A-aI)-\f{A)-f{a)I)B 
/(a) 



hy a function f . 

One can prove an analogous statement in the block-triangular case: 
Proposition 8 Let f he an entire function and let 



Then, 



M 



fiM) 



A B 
C 



a{A)na{C) = 0. 



f{A) f{A)X-Xf{C) 
fiC) 



where X is the (unique) solution to the equation 

AX-XC = B. 
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Proof. Let X be a solution of the Sylvester equation AX — XC = B. Since the spectra of A 
and C are disjoint, this solution is unique [la Section 4.4]. Then, M = T'^ diag(^, C)T where 

/ X 
/ 



Hence f{M) = T ^ diag{ f (A) , f {C))T , which proves the proposition. 



□ 



As an immediate corollary, we obtain an indirect characterization of functions preserving non- 
negativity of block-triangular matrices with two diagonal blocks. 

Corollary 9 An entire function f preserves nonnegativity of block upper-triangular matrices of the 
form 

" A B 



C 



ni xni 



C G IR' 



n2 xn2 



if and only if 

a) f G where N : = maxjni, 71,2}; and 

h) f{A)X-Xf{C) >0 /oreueri/^GlR"!^"!, SelR"^''"^ C € IR^^x"^ such that A, B,C > 0, 
o'{A) n cr(C) = 0, and the (unique) matrix X that satisfies the equation AX — XC = B. 

Proof. For / to preserve nonnegativity of blocks A and C, it has to belong to JF/v (keeping in 
mind Lemma [T]). The remainder of our assertion follows from Proposition [8] and the fact that the 
matrices with nonnegative blocks A, B, C, such that the spectra of A and C are disjoint, are dense 
in the set of all block upper-triangular matrices. □ 

The above proposition, however, does not allow for an explicit formula of the type ([1]) as in 
Theorem [H 

Remark. Note that the results of this section characterize functions preserving nonnegativity of 
the (block) lower-triangular matrices as well. 



6 Preserving nonnegativity of circulant matrices 

A circulant matrix (see, e.g., [7]) A is determined by its first row (ao, • • • , fln-i) as follows: 
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All circulant matrices of size n are polynomials in the basic circulant matrix 
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1 



1 



8 



which imphes in particular that any function f{A) of a circulant matrix is a circulant matrix as 
well. Moreover, the eigenvalues of a circulant matrix are determined by its first row (see ^) by 
the formula 

n— 1 

{ "^^^aj : A; = 0, . . . , n - 1}, where w : =6^^'/". 
Hence the eigenvalues of f{A) are 

n-1 

{fiY^^^'o.,) : fe = 0,...,n-l}. 

j=0 

Thus, the elements (/o, . . . , /n-i) of the first row of f{A) can be read off from its spectrum: 

n— 1 n—1 

= - E ^^'"/(E ^''«^)' ^ = 0, . . . , n - 1. 

fc=o i=o 

This argument proves the following theorem. 

Theorem 10 For an entire function f to preserve nonnegativity of circulant matrices of order n, 
it is necessary and sufficient that for / = 0, . . . , n — 1, 

n—1 n—1 

^ uj-^^fC^ J'^aj) > whenever aj>0, j = 0, . . . , n - 1, where ui = e^'^'/". (3) 

k=0 j=0 

7 Characterization of J^^ for small values of n 

We now focus of the function classes J^n for small values of n. Recall the inclusion J-n+i Q from 
Lemma [H which means that all conditions satisfied by the functions from J^n get inherited by the 
functions from Tn+i- Thus we need to find out precisely how to strengthen the conditions that 
determine Tn to get to the next class J'-'n+i- 

7.1 The case n = 1 

A function / is in J^i if and only if / maps nonnegative reals into themselves. While this statement 
is in a way a characterization in itself, if / is an entire function with finitely many zeros, we can give 
a description of the form that / takes. For such /, the proposition below serves as an alternative 
characterization. 

Proposition 11 A function f having finitely many zeros is in T\ if and only if it has the form 

where the a 's and the (3 's are arbitrary reals, the 7 's are nonnegative, and g is an entire function 
that has no zeros in C and is positive on 1R+. 
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Proof. First note that since / takes real values over the nonnegative reals, all its zeros occur in 
conjugate pairs. Moreover, while the multiplicity of the real negative zeros is not resticted in any 
way, the nonnegative zeros must occur with even multiplicities. This produces exactly the factors 
recorded in ([4]), with nonnegative zeros corresponding to /3 = 0. After factoring out all the linear 
factors, we are left with an entire function — which we call g{z) — that has no zeros, and takes 
only positive values on IR-1-. This gives us the expression ([4]). □ 

Remark. Incidentally, all polynomials / that take only positive values on ]R_|_ are characterized 
by a theorem due to Poincare and Polya (see, e.g., [G', p. 175]): there exists a number G 7Z^ such 
that the polynomial (1 + z)^ f{z) must have positive coefficients. Since we include non-polynomial 
functions into our class Ti, and since we allow functions to have zeros in 1R+, the Poincare-Polya 
characterization is not directly relevant to our setup. 

7.2 The case n = 2 

We just saw that functions in Ti are characterized by one inequality, viz. 

/(x) > Vx>0. (5) 

In this subsection we will see that functions in T2 are characterized by two inequalities, one involving 
a divided difference. We recall two preliminary observations, Lemmas [4] and [5] that were proved in 
Section HI Their specialization to the case n = 2 gives the following corollary. 

Corollary 12 An entire function f belongs to J-2 if and only if it maps positive symmetric matrices 
of order 2 into nonnegative matrices. 

Proof. A strictly positive 2x2 matrix A can be symmetrized by using the transformation 
DAD-^, where D is a diagonal matrix with positive diagonal elements. Thus, Lemmas |4] and [5] 
imply that f{A) is nonnegative for all strictly positive, and hence for all nonnegative matrices A 
of order 2, if and only if f{A) is nonnegative for all symmetric matrices. □ 

Now we are in a position to prove a characterization theorem for the class 
Theorem 13 An entire function f is in T2 if and only if it satisfies the conditions 

f{x + y)- f{x-y) > Vx,y>0, (6) 
(x + y — z)f{x — y) + {z — X + y)f{x + y) > yx>z>0, y>x — z, (7) 

or, equivalently, if f satisfies and the condition 

{x + y)f{x-y) + iy-x)f{x + y)>0 V y > x > 0. (8) 

Proof. If / G then, in particular, / preserves nonnegativity of nonnegative circulant matrices. 
Thus, the conditions ([3]) are necessary for / to belong to J^2- Observe that the condition ([6]) is one 
of the two necessary conditions ([3]) in case n = 2 (taking oq = x and ai = y). Therefore, we need 
to check that the condition ([Tj) is also necessary and that both together are sufficient. Then we 
also need to check that conditions ([6]) and ([7j) are equivalent to conditions ([6]) and ([8]). 
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By Corollary 1121 we can restrict ourselves to the case when ^4 is a positive symmetric matrix, 
i.e., when 

On b 



A 



0.22 



an, 6, > 0. 



Since the value of / at A coincides with the value of its interpolating polynomial of degree 1 with 
nodes of interpolation chosen at the eigenvalues of A p[5i Sections 6.1, 6.2], we get 

f{A) = f[ri]I + f[n,r2]{A-riI), 

where 

,^,='^n+a22^^_^yV{an-a22r+^^ j = l,2. 
So, the off-diagonal entries of f{A) are equal to 

f[ri,r2]b, 

while the diagonal entries are 

f[ri,r2]{ajj - n) + /(n), j = 1,2. 

Writing 

On + 022 



y 



2 ' 

V(aii - 022)2 + 462 



2 

z := min(aii, 022), 

we see that the characterization for J^2 consists precisely of conditions ([6]) and 

It remains to prove that ([6]) and ([7]) are equivalent to ([6]) and By simply taking z = 
in ([7]), we see that dTj) implies ([8]). So let us now assume ([6]) and ([8]). We begin by stating a simple 
auxiliary fact. Taking x = and y > in ([6]) and we get f{y) ± f{—y) > Vy > 0. We 
conclude from this that /(y) > whenever y > — i.e., that / satisfies ([5]). 
First consider y lying in the range x — z < y < x. In this case, we get 

(x + y - z)f(x -y) + {z -x + y)f{x + y) 

= {y-{x- z)){f{x + y)- fix - y)) + 2yf{x - y) > for x > z > 0. 

The nonnegativity of the second term above is a consequence of ([5]), since x — y is nonnegative in 
this case. Now if y > x, then ([6]) and ([8]) simply imply that 

(x + y - z)f{x - y) + (z - X + y)/(x + y) 

= {{x + y)f{x -y) + {y- x)f{x + y)) + z{f{x + y) - /(x - y)) > 
for y > X > z > 0. 

The last two inequalities show that ([6]) and ([8]) imply ([6]) and ([7]). This finishes the proof. □ 
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8 Preserving nonnegative symmetric matrices 



We now focus on the characterization problem for the class of entire functions that preserve non- 
negativity of symmetric matrices. We begin by recalling known facts about functions that preserve 
nonnegative symmetric matrices that are in addition nonnegative definite, i.e., have only nonnega- 
tive eigenvalues. 

8.1 Preserving nonnegative definite nonnegative symmetric matrices 

Interestingly, the condition necessary and sufficient for preserving nonnegative symmetric matrices 
that are nonnegative definite turns out to be exactly the same as the condition for preserving upper- 
(or lower-) triangular nonnegative matrices. 

The characterization of functions that preserve the class of nonnegative definite, entrywise 
nonnegative symmetric matrices is due to Micchelli and Willoughby ^24j. We next state a version 
of their result that is useful for our purposes. 

Result 14 (version of [24, Corollary 3.1]) An entire function f preserves the class of nonneg- 
ative definite, entrywise nonnegative symmetric matrices if and only if all the divided differences 
of f of order up to n are nonnegative over 1R+, i.e., f satisfies (Q]) or, equivalently, all derivatives 
f^^^ of f up to order n—1 are nonnegative on 1R+. 

The proof of Result [T3] in [23] relies on two facts. The first is that f{A) coincides with the 
interpolating polynomial of /, with nodes at the eigenvalues of A, evaluated at A, i.e. that 

f{A) = f[n]I + f[ri,r2]{A - rj) + ■■■ + f[n, . . . r„](A - nl) ■■■{A- rn-iA). (9) 

The second fact is the entrywise nonnegativity of all matrix products 

{A-nI)---{A-rjI), j = l,...n-l, 

which holds under the assumption that the eigenvalues ri, . . . r^ of A are ordered 

ri < r2 < ■ ■ ■ < rn. 

Observe, however, that conditions ([T]) are not sufficient for a function to preserve nonnegativity 
of all nonnegative symmetric matrices. Indeed, let n = 2 and let 

fix) = l+x + ^x'^ -^x^ + ^x\ 

This function satisfies the condition ([T]) with n = 2, but it maps the matrix 

"0 M ' 
M J ' 

which is not nonnegative definite, to a matrix with negative off-diagonal entries when M > is 
chosen to be sufficiently large. In fact, any M > y^3/2 will produce a matrix with negative entries. 

Motivated by Result [lH we would therefore like to find out what conditions are necessary and 
sufficient for a function to preserve nonnegativity of a nonnegative symmetric matrices. We begin, 
in the next subsection, by analyzing even and odd functions. 
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8.2 Even and odd functions preserving nonnegativity of symmetric matrices 



Using the Micchelli-Willoughby result — i.e., Result [U] from the previous section — and an aux- 
iliary result from [12], we shall obtain a characterization of even and odd functions that preserve 
nonnegativity of symmetric matrices. Our proof below will require the notion of a Jacobi matrix 
and that of a symmetric anti-bidiagonal matrix. A Jacobi matrix is a real, nonnegative definite, 
tridiagonal symmetric matrix having positive subdiagonal entries. A matrix A is called a symmetric 
anti-bidiagonal matrix if it has the form 



A 













an 













ai, 



. ,an & K. 



(10) 



and from [12 



a„_2 

O-n Ofi— 1 

We make use of the next two results, from 

Result 15 ([24j) A matrix function f preserves nonnegativity of symmetric nonnegative definite 
matrices of order n if and only if it maps Jacobi matrices of order n into nonnegative matrices 
or, equivalently, if the divided differences of f up to order n satisfy HP for each ordered n-tuple 
Xl < X2 < ■ ■ ■ < Xn of eigenvalues of a Jacobi matrix. 



The above result is not stated in precisely these words in }24j . but it is easily inferred — it lies 
at the heart of the proof of [24] Theorem 2.2]. In addition, we shall also need the following result: 

Result 16 (Corollary 3, |12] ) Let M be a positive real n-tuple. Then, there exists a Jacobi 
matrix that realizes Ai as its spectrum and has a symmetric anti-bidiagonal square root of the 
form U0\) with all aj 's positive. 

We are now in a position obtain a characterization of even and odd matrix functions that are 
of interest to us. 



Theorem 17 An even entire function f{z) =: g{z'^) preserves nonnegativity of symmetric matrices 
of order n if and only if the divided differences of g up to order n are nonnegative on ]R_|_ — i.e., 
if g satisfies (OP. An odd function f [z) =: zh{z'^) preserves nonnegativity of symmetric matrices of 
order n if and only if h satisfies fip. 

Proof. Let / be even. Then, /(z) = g{z'^) for some entire function g. If a matrix A is entrywise 
nonnegative symmetric, then is entrywise nonnegative, symmetric, and nonnegative definite. By 
Result [131 if g satisfies ([1]), then g{A^) is nonnegative. To prove the converse, consider an arbitrary 
n-tuple M of positive numbers. We can think of ^A as being ordered 

M = {Xi, . . . ,Xn), Xi<---<Xn. (11) 

By Result [TUl there exists a nonnegative symmetric anti-bidiagonal matrix A such that A^ is 
a Jacobi matrix with spectrum M. Then, by Result [T5| the divided differences of g must be 
nonnegative when evaluated at the first k points of A^, for each k = 1, . . . ,n. This implies, by the 
standard density reasoning, that all divided differences of g must be nonnegative over ]R-|_. 
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Now let / be odd. Then, f{z) = zh{z'^) for some entire function h. If all the divided differences 
of h up to order n are nonnegative, then by the same argument as above, /i(A^) is nonnegative 
for each symmetric nonnegative matrix A, and multiplication of h^A^) by a nonnegative matrix A 
produces a nonnegative matrix again. To prove the converse, we use induction and a technique 
from [24J. Since / has to preserve nonnegativity of symmetric matrices of order n—1 as well, we 
can assume the nonnegativity of the divided differences of orders k = 1, . . . , n—1. To prove that the 
nth divided difference is nonnegative, let Ai be an arbitrary positive n-tuple (jlip . As above, by 
Result [ini there exists a symmetric anti-bidiagonal matrix A such that is a Jacobi matrix with 
spectrum A4. By [24J, formula ([9]) shows that the {l,n) entry of the function /i(A^) is a positive 
multiple of f[xi, . . . hence the (1, 1) entry of the product h{A'^)A is again a positive multiple 
of f[xi,... ,Xn]- Thus the nth divided difference has to be nonnegative as well, which finishes the 
proof. □ 

This theorem provides a rather natural characterization of even and odd functions that preserve 
nonnegativity of symmetric matrices in terms of their divided differences. However, the "natural" 
idea, that the even and odd parts of any entire function that preserves nonnegativity of symmetric 
functions must be also nonnegativity-preserving, turns out to be wrong. Here is an example that 
illustrates why that may not be the case. 

Example 18 Let 

f{z) ■.=a + I3z - z^ + z^ + 7/, 

where (3 > 1/4, and 0,7 > are chosen to be so large that f{x) > for a// a; E IR and f'{x) > 
for all X € IR+. Then, f preserves nonnegativity of symmetric matrices of order 2, but its odd part 
fodd does not. 

Proof. The function / satisfies conditions ([6]) and ([8]). Indeed, since / > on IR, we have 

{t + s)f{-t) + tf(t + s)>Q Vs,t>0, 

which is equivalent to condition ([8]). Now, the odd part of / is given by 

fodd{z) = l3z-z^ + z^ =: zh{z^). 

Since /3 > 1/4, h{x) > for all x € IR. Since / is monotone increasing on IR+, we have 

/(s + t) - f{-s) > f{s) - f{-s) = 2fodd{s) > V s, t > 0, 

which yields condition ([6]). Thus, by Theorem 1131 / preserves nonnegativity of symmetric matrices 
of order 2. However, 

h'{x) = 2x - 1< for X < 1/2. 

Therefore, by Theorem 1171 fodd does not preserve nonnegativity of symmetric functions of order 2. 
□ 

We conclude this section with a simple observation about even and odd parts of a nonnegativity- 
preserving function. 

Proposition 19 // an entire function f preserves nonnegativity of symmetric functions of order 
n, then its odd and even parts fodd o^'^d feven preserve nonnegativity of matrices of order [n/2j . 



14 



Proof. For n even, consider matrices of the form 

A : 

and for n odd, matrices of the form 

^ = diag( ^ ^ ,0), 
where B is an [n/2j x [n/2j symmetric nonnegative matrix. Since 

for n even, 



B 
B 



m = 

f{A) = diag( 



feven{B) fodd{B) 
fodd{B) feven{B) 



, 0) for n odd. 



feven{B) fodd{B) 
fodd{B) feven{B) 

the see that feven and fodd must preserve nonnegativity of symmetric functions of order [n/2j. □ 



8.3 Other necessary conditions 



Results from |12j allow us to derive an additional set of necessary conditions. The motivation 
behind these conditions is as follows. We believe that the power of Results [15] and [16] — or rather, 
the methods behind those results — have not been exhausted by Theorem 1 171 Our next theorem is 
presented as an illustration of this viewpoint. On comparison with Theorem 1131 we find that the 
conditions derived in our next theorem constitute a complete characterization for the functions of 
interest in the n = 2 case. To derive these new necessary conditions, we will need the following two 
results. 

Result 20 (Theorem 1, [12]) A real n-tuple A can be realized as the spectrum of a symmetric 
anti-bidiagonal matrix U0\) with all aj 's positive if and only if A = (Ai, . . . , A„) satisfies 

Ai>-A2>A3>--->(-l)""^A„>0. 

Lemma 21 Let A be a symmetric anti-bidiagonal matrix of order n, and let Afj denote the {i,j) 
entry of A^. Then 

a) The {i,j) entry of A'^'^"^ is zero whenever 2 <i+j < {n — q + 1), q > 1. 

b) The {i,j) entry of A?"^ is zero whenever 1 + (7<J — i<n — 1, g>l. 

c) Adopting the notation in U0\) for the entries of A, 

'^''"-^ - OnOn-l . . . an-2q+2, 1 < Q < [(n + l)/2j , 

anttn-l ■ ■ ■ an-2q+l, 1 < 9 < [1^/2]. 



^l,n-q+l 
^1,1+9 



(12) 
(13) 
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Proof. We proceed by induction on q. Note that (a), (b) and (c) are obvious when q = 1- Let 
us now assume that (a) and (b) are true for some q < n — 3. Note that since A is anti-bidiagonal, 

^ii^^ = A>-j+i^"-i+ij + ^i,n-i+2^"-i+2j- (14) 
However, if i + j < (n — (g + 1) + 1), then 

{n - j + 2) - i > [n - j + 1) - i > q + I. 

Applying our inductive hypothesis on (b), we conclude from the above inequalities that the right- 
hand side of (fTll) reduces to zero when i + j < (n — (g + 1) + 1). Thus, (a) is established for 
q + l. 

We establish (b) for g + 1 in a similar fashion. We note that 

A2q+2 _ .2q+l . _L/l2g+l . 

When j — i > 1 + (g + 1), then 

i + {n - j + 1) < i + {n - j + 2) < n - (q + 1) + 1. 

Since we just established (a) for q + l, the above inequalities tell us that the right-hand side of ([15]) 
reduces to zero when j — i > 1 + (g + 1). Thus, (b) too is established for g + 1. By induction, (a) 
and (b) are true for all relevant g. 

Part (c) now follows easily by substituting i = 1 and j = n — q into equation (I14p to carry out 
the inductive step for (fT2]) , and by substituting i = 1 and j = g + 2 into equation (fTSl) to carry out 
the inductive step for (fT3|) . □ 

We can now present the aforementioned necessary conditions. 

Theorem 22 // an entire function f preserves nonnegativity of symmetric matrices of order n, 
n>2, then, for each ordered n-tuple {xi, . . . , Xn) where 

Xi > -X2 > X3> ■■■ > (-l)""^x„ > 0, (16) 

/ must satisfy 

f[xi,...,Xn]>0, (17) 

and for each k = 1, . . . ,n, f must satisfy 

f[xi, . . .,Xk-l,Xk+l, ...,Xn\-{ ^Xj)f[xi, . . . > 0. (18) 

Proof. We choose an n-tuple (xi, . . . ,x„) that satisfies (fT6]) . By Result \20\ there is a symmetric 
anti-bidiagonal matrix of the form (jlOp . with all Oj's positive, whose spectrum is (xi, . . . ,x„). Let 
us express f{A) using the formula Q, with the substitutions rj = Xj, j = 1, . . . ,n. Then, in view 
of Lemma [2T| the (1, [n/2j + 1) entry of f{A) is a„a„_i . . . a2/[xi, . . . , Xn]- Since / preserves non- 
negativity, and all the Oj's are positive, f[xi, . . . ,Xn] has to be nonnegative. This establishes (fT7|) . 
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To demonstrate (jlSp . we look at the entries of f{A) that are adjacent to the (1, \ n/2\ + 1) entry 
that was considered above. Let us fix a A; = 1, . . . ,n. This time, however, in using formula ([9]) to 
express f{A), we make the following substitutions 

{Xj if j < k, 
Xj^i ii k < j < n, 
Xk if j = n. 

Our analysis splits into two cases. 

Case 1. n is odd: In this case, let us look at the (1, \ji/2\ + 2) entry of f{A). By Lemma [2H 
and the fact that n is odd, the only power of A that contributes to this entry is A^~'^. Consequently 

/(^)l,Ln/2j+2 = {f[xi, ■ ■ .,Xk-l,Xk+l, ...,Xn\-{ ^Xj)/[xi, . . • , a^n]}^" J^/2j+2 

= anan-i ■ ■ -a^iflxi, . . . ] - ( '^Xj)f[xi, ■ ■ ■,Xn\]. 

Since / preserves nonnegativity, (jlSp follows from the above equalities. 

Case 2. n is even: In this case, we focus on the (1, \ n/2\) entry of f{A). We recover (llSp by 
arguing exactly as above. 

In either case, (jlSp is established, which concludes our proof. □ 

We conclude this section by showing that a subset of the necessary conditions derived above 
are in fact sufficient to characterize those entire functions that preserve nonnegativity of 2 x 2 
symmetric matrices. Specifically, we show that 

f[xi^X2\ > and 
f{x2)—X2f[xi,X2] > for all xi > — X2 > 0, 

imply the conditions ^ and ([8]). This is achieved simply by taking some y > x > 0, making the 
substitutions xi = y + x and X2 = x — y, and then invoking continuity to obtain ([6|) and ([8]) for all 
y > X > 0. 

9 Open problems and further ideas 

We conclude this paper by listing some ideas that we did not pursue, which however may lead to 
further progress. 

One can consider matrices that preserve nonnegativity of other classes of structured matrices, 
such as Toeplitz or Hankel. However, since these classes are not invariant under the action of 
an arbitrary matrix function, their matrix functions can be quite difficult to analyze. Also, the 
eigenstructure of some structured matrices is rather involved, which could be an additional obstacle. 

Theorem 1.3 of |35j gives an interesting formula for f{A) when / is a polynomial, which therefore 
must also be true for entire functions. Precisely, if ^ is a matrix with minimal polynomial po and 
C is the companion matrix of then 

n 

f{A)=Y^f{C)j,A^-\ 
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In particular, f{A) is nonnegative whenever the first eolumn of /(C) is nonnegative. It would be 
worthwhile to find out what functions have this property. 

Note that the set J^n contains positive constants and is closed under addition, multiplication, 
and composition. We are not aware of any work on systems of entire functions (or even polynomials) 
that satisfy this property. Perhaps one could describe a minimal set of generators (with respect to 
these three operations) that generate such a system. 

For example, in the case n = 1, the generators are positive constants, the function pi{x) = x 
plus all quadrics of the form {x — of', a > 0. Incidentally, the set of polynomials with nonnegative 
coefficients is generated by positive constants and p\{x) = x. We do not have a characterization of 
generators for n > 2. 

In particular, J^n is a semigroup with respect to any of these operations, so some general results 
on semigroups may prove to be useful in our setting. Also note that the set of nonnegative matrices 
of order n, on which Tn acts, is also a semigroup (closed under addition and multiplications), which 
could also be of potential use. 

Finally, both J^^i and the set of nonnegative matrices of order n are also cones, so the problem 
might also have a cone theoretic form. If we consider polynomials instead of entire functions, we 
can further restrict ourselves to polynomials of degree bounded by a fixed positive integer. Then, 
we will obtain a proper cone, whose extreme directions may be of interest. The general problem 
then can also be looked upon in an appropriate similar setting. 
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